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ABSTRACT 

Although an understanding of radiometric dating is central to the preparation of every geologist, many students 
struggle with the concepts and mathematics of radioactive decay. Physical demonstrations and hands-on experiments 
can be used to good effect in addressing this teaching conundrum. Water, heat, and electrons all move or flow in 
response to generalized forces (gradients in pressure, temperature, and electrical potential) that may change because of 
the flow. Changes due to these flows are easy to monitor over time during simple experiments in the classroom. Some 
of these experiments can be modeled as exponential decay, analogous to the mathematics of radioactive decay, and can 
be used to help students visualize and understand exponential change. Other, similar experiments produce decay or 
change that is not exponential. By having classes, in small groups, conduct several experiments involving flows, a 
learning synergy can be encouraged in which the physical and mathematical similarities of flow processes are 
emphasized. For the best results, students should be asked to analyze the experimental data, using graphs and algebra 
or calculus as appropriate to the class, to determine the nature of the decay process and to make predictions, either 
forward or backward in time as would be done for radiometric dating. Basic quantitative skills are strengthened or 
developed as part of these activities. Encountering a number of important geologic processes in the same mathematical 


INTRODUCTION 

Radiometric dating provides essential data for much 
of our geologic understanding of earth history and, 
therefore, is an important topic for many courses in 
geology. Because calculus is typically used in developing 
the mathematics of radioactive decay, students without 
calculus are likely to have difficulty in grasping some of 
the principles of radiometric dating and in applying the 
relevant equations. The significance of radioactivity, 
exponential decay, and radiometric dating to geologists, 
chemists, and physicists and the challenges of teaching 
about these topics have led to numerous published 
pedagogic strategies. These include stochastic experiments 
with coins (Kowalski, 1981; Tyburczy, 2000; Wenner, 
2008a), dice (Celnikier, 1980; Priest and Poth, 1983; 
Schultz, 1997; Benimoff, 1999), poker chips (McGeachy, 
1988), M&M's (AAAS Science Netlinks, 2001a; Gardner et 
al., 2005; Wenner, 2008b), or popcorn (Wenner, 2008c), 
measurements of actual radioactive decay (Jones, 1957; 
Whyte and Taylor, 1962; Supon and Kraushaar, 1983), 
analog experiments (Bohn and Nadig, 1938; Knauss, 1954; 
Smithson and Pinkston, 1960; Wunderlich and Peastrel, 
1978; Wise, 1990; AAAS Science Netlinks, 2001b; 
Greenslade, 2002; Fairman et ah, 2003; Sunderman, 2007), 
and general mathematical approaches or advice 
(Guenther, 1958; Vacher, 2000; Shea, 2001; Huestis, 2002). 
I have tried a number of classroom activities to help 
students understand the mathematics of radioactive 
decay, and describe here some of the more successful 
experiments that I find useful for geology classes. I 
believe that having classes do several of the experiments, 
simultaneously in small groups, demonstrates the 
mathematical connectedness of a variety of important 
geologic processes. Furthermore, the data collection and 
analysis skills learned by studying any one of the 
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processes help build a quantitative literacy that can be 
transferred and applied to the other processes. 

A SIPHON EXPERIMENT 

My most-used decay experiments involve flowing 
water, and can be done with a "U-tube" or a siphon. 
Geologists think a lot about flowing water. Streams, 
erosion, groundwater, and flooding are topics in 
introductory and advanced geology courses, and some of 
the principles of flowing water are familiar to even 
beginning students. I find that most geology students are 
able to make good predictions when presented with the 
experimental design shown in Figure 1. The "U-tube" 
apparatus consists of two vertical plastic tubes (44.5 mm 
i.d.) connected at their bases by a copper pipe (12.7 mm 
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i.d.) and valve. The valve is designed so that it can be 
opened quickly to a preset p>osition against a setscrew (see 
inset in Figure 1). What will happen if the tube on the left 
is filled with water to the height Iil, the tube on the right is 
filled with water to the height hR, and the valve is 
opened? All students predict correctly that water will 
flow from the mostly full tube to the mostly empity tube 
through the connecting pipe until the water level in the 
two tubes is equal. And if asked, many students will 
predict that the flow will slow down as the difference in 
water heights (Iil-Iir) decreases. Their intuition, based on 
previous experiences, tells them (or will let them be easily 
convinced) that the flow rate is prop>ortional to the 
difference in heights, analogous to hydrologic head 
ill = hi-huj: 

Water Flow = -krH (1) 

where ki is a constant. (The minus sign is used here to 
anticipate other flow equations where a coordinate system 
is used. The subscript 1 in the constant ki is used to 
differentiate ki from the other constants ki used in this 
paper.) 

The plastic tubes are cylinders, so the volume of water 
flow is proportional to the change in head (dH). Replacing 
"Water Flow" by dH/dt, equation (1) becomes 

dH/dt = -k 2 -H (2) 

which is an exponential decay equation. The constant k 2 
in (2) is related to the constant kj in (1), but incorporates 
the volume to height conversion. Because of its 
mathematical similarity, the change in H with time is a 
visual analog for the change in parent (or daughter) 
isotopes during radioactive decay. During the experiment, 
the height of the water in either tube can be seen to change 
quickly at first, then more and more slowly as the water 
heights in the two tubes approach one another. A little 
food coloring makes the water easier to see. Note that 
dhi.= -dhR, so that 

dH/dt = 2 dlu/dt = -2 dh R /dt = -k 2 -H (3) 

I have used this experiment both as a demonstration 
in front of a large audience and as an experiment for a 
small group in a lab setting. By adjusting the valve 
opening, the complete experiment can be done in less than 
two minutes, or over a longer time p>eriod to allow for 
careful data collection of water heights as a function of 
time. Although the U-tube apparatus in Figure 1 is not 
difficult to build, the experiment can be done with 
minimal preparation by using two large graduated 
cylinders joined by a clear pilastic sipihon tube. I used the 
siphon setup for several years before building the Figure 1 
U-tube device for convenience, so I will always think of 
this as a siphon experiment. And after discovering the 
papers by Knauss (1954) and Smithson and Pinkston 
(1960), 1 have started using a siphon setup again for 
student data collection, as 1 will explain shortly. 

How much mathematical detail I use and whether I 
have students collect data on the change in water height 
with time depends on my class and goals. In most cases, I 
present the U-tube apparatus with one tube full of colored 
water, discuss with the class their predictions about the 


consequences of opening the valve, and then proceed with 
a short experiment. I then develop the mathematics of 
exponential decay, beginning with an algebraic derivation 
of the decay equations (given later in this paper) that 
doesn't involve calculus. The recursive character of the 
algebraic approach is well suited for adapting to a 
spreadsheet model and graph, with each line in the 
spreadsheet representing one unit of time (e.g. one half- 
life). Students in some settings can be asked to do this as 
an exercise, during or between classes, or by the instructor 
using a compu ter and projector. After obtaining the decay 
equation (28), I then show the calculus derivation for those 
students who have had calculus. With an exponential 
decay equation in hand, I may revisit and p>erhap>s rerun 
the U-tube (or siphon) experiment. The concept of a half- 
life is very useful here and transfers directly to the visual 
demonstration involving water heights. 

1 typically have my students collect their own data on 
water height as a function of time for the U-tube or siphon 
setup), either as a whole class experiment or as a small- 
group) lab assignment. My preferred approach is to have 
each small group* collect data on a different process or 
experimental setup), such as various water flow designs, 
cooling objects, or capacitor discharge. I describe some of 
these other experiments later in this papier. Most are 
exponential decay experiments, but some are not and 
serve to demonstrate the power of data analysis. The 
assignment includes a lab report that must include (1) a 
presentation of the data in graphical as well as tabular 
form, (2) a calculation of half-life and A (decay constant) 
for their experiment, and (3) an evaluation of uncertainty 
in their measurements and calculations. 

Data from a typucal 10-minute-long U-tube 
experiment using the Figure 1 apparatus are shown in 
Figure 2A. Shown also in Figure 2A is a solid line 
grapihing exponential decay for equation (28) using a 
value of A- taken from the first half-life observed (the time 
required to reduce H to one-half its starting value). 
Similar results are obtained for a siphon experiment using 
a 6 mm i.d. plastic siphon tube with a clamp to adjust flow 
rate. Note that the time required to equalize the water 
levels in the two tubes is less than predicted by the first 
half-life. Interestingly, the observed half-lives get shorter 
and shorter during the experiment. Although some 
choices of A fit the data better than the value from the first 
half-life, there is no single value of A that will describe an 
exponential decay that matches all the data. This 
discrepancy can be ignored by using an average value of 
A, or it can used as a point of discussion and perhaps 
further experiments. A plausible explanation for the 
changing half-life length is that the resistance to flow of 
the water is greater when the water velocity is higher. 
This hypothesis can be tested by reducing the U-tube flow 
rate with the connecting valve (or by pinching the siphon 
tube with a clamp). Results of longer experiments are 
shown in Figure 2B. To aid in comparison, the data have 
been normalized by dividing the H values by Ho and the 
time values by the observed first half-life. Results from 
slower experiments do indeed more closely approximate 
the exponential decay equation, supporting the notion 
that resistance due to turbulence in the valve is 
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responsible for the observed decrease in half-life as H 
decreases. 

Smithson and Pinkston (1960), following Knauss 
(1954), present a hydrodynamic design that yields good 
exponential decay data rising a tall vertical tube and a 1.25 
-meter-long capillary tube drain, making the flow slow 
enough to be viscosity-controlled and follow Poiseuille's 
Law (see also Skinner, 1971). Their results encouraged me 
to try a 1.25-meter-long capillarv-sized (2 mm i.d.) plastic 
siphon tube, instead of the larger (6 mm i.d.) tubing I had 
been using. Data from a siphon experiment with the 2 
mm i.d. siphon are shown on a semi-log graph in Figure 
2C, along with some of the U-tube data of Figure 2B. If In 
(WHo) is plotted vs. time t, exponential decay data should 
fall along a straight line with slope - X (see equation 27). 
The siphon data fall nicely along the exponential decay 
line in Figure 2C, showing the advantage of the small 
tubing size and absence of a constriction. A drawback is 
the length of the experiment, but this can be mitigated by 
using smaller i.d. vertical tubes or by collecting only part 
of the data, which would still permit an evaluation of X. 
Because of the variability of the siphon start, results are 
better if you get the siphon flowing before recording the 
starting value of H (=Ho) for your data set. 

A related experimental option also may be available if 
someone in your department studies groundwater. A 
common tool for laboratory measurement of hydraulic 
conductivity is a permeameter, which may be used in a 
constant-head in or a falling-head geometry (Fetter, 2001). 
A falling-head experiment is similar to the Smithson and 
Pinkston (1960) experiment, but with a soil sample 
replacing the capillary tube in the role of slowing the flow 
and keeping it in the lamellar flow regime. Some effort 
and experimentation is needed to find a soil sample that 
has a low enough permeability to allow precise 
measurement of head height as a function of time. 


FIGURE 2. Data from U-tube and siphon experiments. 
(A) Data from a 10-minute-Iong U-tube experiment are 
shown in terms of the height difference FI and time t. 
Ideal exponential decay is shown as a solid curve using 
the experimental time for H to fall to half its starting 
value (144 seconds for this experiment) to define the half- 
life length. Note that the second and subsequent half- 
lives of the siphon experiment are shorter than the first 
half-life. (B) Data from three U-tube experiments are 
shown in terms of the fraction H/Ho of the original height 
difference Ho and time t given as "half-lives", using the 
experimental time for H/Ho to fall to 0.5 to define the half 
-life length for each experiment. Ideal exponential decay 
is shown as a solid curve. Longer siphon experiments 
with slower flow rates more nearly approximate ideal 
exponential decay. (C) Data from a siphon experiment 
using a 2 mm i.d. siphon are shown in terms of the natu¬ 
ral log of the fraction H/Ho of the original height differ¬ 
ence Ho and time t given as "half-lives", using the first 
experimental half-life to normalize the data. Ideal expo¬ 
nential decay is shown as a solid straight line. Also 
shown for comparison are U-tube experimental data from 
(A). The siphon data follow an exponential decay trend 
very closely when a "capillary-sized" siphon is used. 



FIGURE 3. Drawings of water clock designs discussed in 
the text. (A) A cylinder with a hole near the bottom. (B) A 
cone with a hole at the bottom and water draining into a 
cylindrical collection vessel. (C) A bucket-shaped water 
clock with a bottom that is one-half the diameter of the 
top. (D) Height of the water h in a cylindrical tank that is 
leaking from a hole near the bottom (A) shown as a func¬ 
tion of time t. The height h is given as a fraction h/ho of 
the starting height ho and the time t is given in units of 
half-lives, using the time for the water level to fall to one- 
half the original height to define the half-life. Ideal expo¬ 
nential decay is shown as a dashed curve. Neither the 
leaking tank nor exponential decay exhibits linear varia¬ 
tion of water height with time, especially after the first 
half-life. 
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However, a falling-head permeameter will also yield good 
exponential decay data and its use makes a very direct 
connection for students of the utility of the mathematics 
for hydrologic studies. 

WATER CLOCKS 

In many respects, the U-tube and siphon experimental 
apparatus described here can be considered water clocks: 
flowing water leads to a change in the water level that 
varies systematically with time. However, the exponential 
decay character of the siphon water clock would have 
been inconvenient for those wishing to track time in the 
ancient world. Water clocks were used lav Egyptians at 
least 3500 years ago and bv other ancient cultures 
including those in Greece (where they were called 
clepsydras, which means "water thief"), the Middle East, 
and China (Sloley, 1931; Barnett, 1998). Many students 
will have learned about water clocks in other contexts and 
may have questions concerning the relationship of the 
siphon experiment to ancient water clocks. Indeed, some 
students may have previously constructed a water clock, 
as there are many sources describing water clock 
experiments for elementary and middle school students 
(e.g. Zubrowski, 1988; Grossman et al., 2000; National 
Geographic Kids, 2008). 1 recommend capitalizing on this 
previous connection and including water clocks as one of 
the options for small group experiments. For some 
classes, it may even be a good strategy to begin with 
traditional water clocks and to discover or confront the 
behavior of the "siphon clock." For these reasons, I 
describe here the interesting quantitative features of 
several simple water clocks. 

Ancient water clocks came in many designs, but most 
were based on water dripping or squirting out of a hole in 
the bottom or side of a vessel, some with the flow 
collected in another vessel (see Figure 3 for some water- 
clock designs). 

Time was measured by the fall in the water level of 
the draining vessel or by the rise in the water level of the 
collection vessel. Although the rate of flow out of the 
ancient water docks does depend on the level of water in 
the draining vessel, it turns out that drainage does not 
lead to an exponential decay of the water level. 

Water draining from a cylindrical tank with a hole in 
the bottom follows Torricelli's Law (Mills, 1982; Groetsch, 
1993), named for the mathematician Evangelista Torricelli 
who lived in Italy in the 17 lh century. The kinetic energy 
of the squirting water changes the flow dynamics and 
must be accounted for in the mathematical description. 
Torricelli's Law states that the flow will be proportional to 
the square root of the water level: 

(^) = <W)G <4) 

where V is the volume and h is the height of water in the 
tank, t is time, and ks is a constant. Equation (4) can be 
derived from Bernoulli's equation for energy balance in an 
incompressible liquid (Ansaldo, 1982; Fetter, 2001, pi.115). 
Because V=h(jtr 2 ) for a cylinder of radius r, we can replace 
V with h in equation (4) and constant ki with constant 
k4=k-/(7rr 2 ) to yield: 


fdh' 

1 dt j 


( - k ^ w/r ( i \,w 

(- r)3/h = (-k 4 Wh- 

71 • r“ 


( 5 ) 


Integration of (5) yields an expression for the water level 
in the cylindrical tank as a function of time: 

h J2zV0f (6) 

4 

assuming that h=l at t=0. The interested reader can verify 
this result by differentiating (6) with respect to t and 
comparing the result with (5), using (6) to get the value of 
\'h . A graph of equation (6) is shown in Figure 3D, along 
with an exponential decay graph for a process with the 
same first half-life. See also Farmer and Gass (1992) and 
Driver (1998). 

It is apparent from Figure 3D that the variation of 
water level in a cylindrical water clock is not much 
improvement (if any!) on the siphon water clock, if a 
linear variation with time is sought. One solution to this 
problem was to keepi the water level constant in the 
dripping vessel by constantly supplying water, adding an 
overflow, and measuring time from the level of water in a 
cylindrical collection vessel. A more common and simpler 
solution was to build water clocks with sloping sides, in 
effect replacing the cylinder with a cone. This makes sense 
because more water must flow from the vessel when the 
cone is full and the flow rate is highest. But does the 
water level in a cone-shaped water dock vary linearly 
with time? Does the cone angle matter? 

For a cone, the volume V is related to the height h and 
the half-angle <(> by 


V - 


| (tan. ^6) h 3 


(7) 


Differentiating (7 with respect to h and using (4) we find: 

(-k-J Vh 


fdh' 


' dh ' 

( dv' 

V dt 7 


v dV ) 

i dt J 


7i (tan 2 <j>)h 2 


( 8 ) 


Solving (8) for a cone-shaped, dripping water clock 
with an initial water depth of one (h=l for t=0) gives the 
expression for the water height h as a function of time t: 


2 5 k?, t 

\ 2 n (tan ~ 4») / 


(9) 


This equation is plotted in Figure 4A for two cones of 
different half-angles <)>, along with a graph of h for a 
cylinder, choosing values of ka so that each vessel drains 
completely in the same time period. Cones with larger 
cone angles hold more water for the same height, so the 
flow rate and therefore the value of k 3 must be greater to 
drain in the same time as a cone with a smaller cone angle. 
Note, however, that the normalized result is independent 
of the cone half-angle! Interestingly, the variation of water 
level h with time for the cone water dock is not much 
improvement (if any!) on the cylindrical water dock or the 
siphon dock. 

The real value of the cone shape for a water dock 
derives from the fact that the flow rate (dV/dt) is nearly 
independent of time. Substituting for h using (9) in 
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FIGURE 4. Results of water clock models. (A) Height of the water h in a cone-shaped water clock (Figure 3B) that is leaking 
from a hole at the bottom of the cone shown as a function of time t. The height h is given as a fraction li/ho of the starting height 
ho and the time t is given as a fraction t/tF of the total time tF needed to drain the water clock. Calculated results for cone-shaped 
water clocks with very different cone angles (10° and 80°) are identical when normalized in this way. Also shown for comparison 
is the time variation of the water height in a cylindrical water clock (Figure 3A). (B) Volume of the water V of the water in a cone- 
shaped water clock (Figure 3B) that is dripping from a hole at the bottom of the cone shown as a function of time t. The volume V 
is given as a fraction V/Vo of the starting volume Vo and the time t is given as a fraction t/tF of the total time tp needed to drain the 
water clock. Also shown for comparison is the time variation of the volume of water in a cylindrical water clock (Figure 3A). The 
volume of water released from the cone-shaped water clock exhibits a nearly linear variation with time over a large fraction of the 
time needed to drain the clock. Therefore, the height of water in a cylindrical collection vessel would rise uniformly with time. 

(C) Height of the water h in a bucket-shaped water clock (Figure 3C) that is leaking from a hole at the bottom of the bucket is 
shown as a function of time t. The height h is given as a fraction h/ho of the starting height ho and the time t is given as a fraction 
t/tF of the total time tp needed to drain the water clock. Calculated results for bucket-shaped water clocks with bottoms that are 
different fractions of the height of a cone are shown. The height of the water in a bucket-shaped water clock that has a bottom at 
0.5 exhibits a nearly linear variation with time over a large fraction of the time needed to drain the clock. Also shown for compari¬ 
son is the time variation of the water height in a cylindrical water clock (Figure 3A) and a cone-shaped water clock (Figure 3B). 

(D) Volume of the water V in a bucket-shaped water clock (Figure 3C) that is leaking from a hole at the bottom of the bucket 
shown as a function of time t. The volume V is given as a fraction V/Vo of the starting volume Vo and the time t is given as a frac¬ 
tion t/ti of the total time tr needed to drain the water clock. Calculated results for bucket-shaped water clocks with bottoms that 
are different fractions of the height of a cone are shown. Also shown for comparison are the time variations of the volume of wa¬ 
ter in a cylindrical water clock (Figure 3A) and a cone-shaped water clock (Figure 3B). Of these water clocks, the rate of water re¬ 
lease from the cone-shaped clock is the most uniform with time. 


equation (7), keeping the initial water depth of one (h=l 
for t=0), V is given by: 


V = 


(2 Jt - 5 k 3 t cot up) 1 


1 tan 2 ([) 


6 (2ji) 


i / = 


( 10 ) 


Equation (10) is plotted in Figure 4B, where the 
volume fraction of water V/Vo that has dripped from a 
conical water clock is shown as a function of time (Vo is 
the initial volume of water in the cone). The volume of 


water dripped from a cylindrical water clock is shown for 
comparison. It is evident from Figure 413 that the rate that 
water that drips from a cone-shaped vessel is very nearly 
constant with time, so the rise of the level of water in a 
cylindrical collection vessel is steady and makes a 
reasonable clock. 

In spite of this good result for two-vessel, cone-shaped 
water clocks, many examples ol one-vessel water clocks 
are known. These have sloping sides, but are only the top 
part of a cone with a flat bottom and the shape of a 


200 


Journal of Geoscience Education, v. 57, n. 3, May, 2009, p. 196-205 















modern bucket (see Figure 3C). I was unable to obtain an 
analytical solution for the water level h as a function of 
time t for a "bucket" water clock, but numerical solutions 
can be obtained with a spreadsheet using equations (4) 
and (6). Results for h and V over time from the numerical 
tests of "bucket" water clocks are shown in Figures 4C 
and 4D. Data are shown for "buckets" with bottoms at h 
= 0.2, 0.5, and 0.8 relative to a full cone height of 1. 
Buckets with a bottom at h = 0.2 are nearly a full cone and 
give results similar to cone-shaped water clocks. Buckets 
with a bottom at h = 0.8 are more like a cylinder and give 
results similar to a cylindrical water clock. Interestingly, 
buckets with a bottom at h = 0.5 balance the features of a 
cone and a cylinder to create a water dock that has a 
nearly linear variation of h with time. The h = 0.5 bucket 
clock is very close to the shape of a water clock in the 
Cairo Museum that dates to the reign of King Amenhotep 
III (1415-1380 BC) (Science Museum, 2008). Was this the 
result of a lucky guess? Good mathematical analysis? Or 
many experiments? 

Water clocks can be made from vessels of most any 
shape, with a decay/draining behavior that reflects the 
shape. You can provide your students with vessels to use, 
or you can let them find their own. Funnels can be used 
for cone-shaped water docks, but the size of the hole in 
the bottom will need to be reduced by a plug, which can 
be a cork or even chewing gum. Vessels for bucket water 
clocks with a 0.5 bottom height are not easy to find; most 
common containers have a bottom height over 0.6. A 
plastic funnel can be used as a "bucket" by plugging the 
bottom and drilling a drain hole at the desired bucket 
bottom height. When you find a vessel you want to use, 
start with a small drainage hole, enlarging it if you need to 
increase the drainage rate. Depending on the container, 
you may be able to have multiple drain holes with plugs 
in, or duct tape over, the ones not in use. An excellent 
discussion of various water clock designs and results of 
some experiments to test "linear outflow clepsydra" can 



be found in Mills (1982). 

COOLING COFFEE 

Heat transfer is another process that can be described 
by the mathematics of exponential decay. The flow of 
heat by conduction between an object (such as a cup of 
coffee) and its surroundings is proportional to the 
difference in temperature (T-Ts) 

Heat Flow = -ks/T-Ts) (11) 

where T is the temperature of the object, Ts is the 
temperature of the surroundings, and ks is a constant 
(Fourier's Law). If the object and its surroundings are 
"well stirred" or the heat flow within the object and 
within the surroundings is rapid relative to the heat 
transfer between them, the gain or loss of heat leads to a 
change in temperature of the object that depends on the 
heat capacity of the object. The resulting change in 
temperature dT/dt, is directly related to the heat flow (11) 
so that: 

dT/dt = -k„-(T-T s ) (12) 

which is an exponential decay equation. The constant kf, 
in (12) is related to the constant ks in (11) with heat 
capacity and surface area terms added. Because many 
students will have had exp>erience with cooling coffee or 
tea or cocoa, measurement of temperature as a function of 
time for coffee is another good small group experiment to 
include in your decay mix. 

Many physicists, mathematicians, and engineers have 
used cooling coffee in their teaching, and a detailed 
analysis can be quite complicated (e.g. Walker, 1977; 
Dennis, 1980; Rees, 1988). Of particular interest has been 
the question of when to add cream if one wishes to speed 
the cooling of the coffee to drinking temperature (e.g. 
Greenslade, 1994; Smith, 2008). I have found cooling 
coffee to work well as an exponential decay experiment, in 
spite of these complications. We fill a paper cup with hot 
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FIGURE 5. (A) Experimental data for cooling coffee in 12-ounce paper cups, uncovered and covered. The logarithm 
of the fractional difference between the temperature of the coffee T and the temperature of the surroundings Ts is 
shown in terms of time given in cooling half-lives. To is the measured temperature of the coffee when the cooling 
started (uncovered 82°C, covered 92°C). The data for the two experiments were normalized based on their own tem¬ 
peratures and first half-lives (uncovered 28 minutes, covered 54 minutes). The covered cup data more closely ap¬ 
proach exponential cooling, which is shown by the straight line. (B) Circuit diagram for a simple capacitor discharge 
experiment. Two D batteries in series are used to charge the capacitor by closing the bottom switch. The capacitor is 
discharged by opening the lower switch and closing the upper switch. As noted in the text, the voltage of the capaci¬ 
tor will decrease exponentially with time and can be monitored with the voltmeter. 
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water (or coffee or tea), cover it with a lid (sealing the 
original holes if any), then insert a thermometer through a 
new hole in the center of the lid using a rubber stopper on 
the thermometer to preset the insertion length. The lid not 
only holds the thermometer, but it also prevents 
evaporation of the hot water and the attendant heat loss, 
which is significant. Figure 5A shows normalized cooling 
curves for covered and uncovered paper cups. The 
uncovered cup loses heat more rapidly at higher 
temperatures, principally because of evaporation. 
Covered or not, the paper cups lose heat less rapidly at 
lower temperatures than predicted by exponential decay 
based on the first half-life. This may be due to less 
vigorous convection in the colder water, which could 
make the water temperature less homogeneous. The 
cooling half-life of a 12-ounce (400 ml when completely 
full), covered paper cup is over 50 minutes, so data 
gathering limited to an hour will yield a result that is close 
to exponential decay. Ceramic mugs have cooling half- 
lives similar to paper cups of the same volume. 

An alternative simple cooling experiment that yields 
more nearly ideal exponential decay data over several half 
-lives was described by Dewdney (1959). He used metal 
cylinders with a hole drilled in the center of the end that 
was large enough to insert a glass thermometer to the 
middle of the cylinder. I have tried this with a 38.1 mm 
o.d. aluminum cylinder, 76.2 mm long with a 6.6 mm 
diameter axial hole for our red-liquid, glass thermometers. 
This cylinder has a cooling half-life of 20 minutes and is 
very easy to use. Because the thermal conductivity of 
aluminum is high relative to the rate of heat loss to the 
surrounding air, the "well-stirred" boundary conditions 
are satisfied until the temperature approaches room 
temperature. 

CAPACITOR DISCHARGE 

Perhaps the most ideal analog experiment for 
exponential decay is to measure the charge on a capacitor 
as it discharges through a resistor. This experiment does 
not provide the visual display that accompanies the 
siphon experiment (unless you add a light bulb to the 
circuit), or the familiarity of cooling coffee. Nevertheless, 
the experiment is easy to set up and you can design a 
layout that meets your half-life needs. The current I 
(ampere = coulombs/s) in a simple dc-circuit consisting of 
a capacitor and a resistor (see Figure 5B) obeys Ohm's law 

I = V / R (13) 


where R is the resistance (ohms) and V is the voltage 
(volts). The voltage V of a charged capacitor is given by 
the charge Q (coulombs) divided by the capacitance C 
(farad = coulomb/volt), and the current I is equal to the 
change of charge with time dQ/dt. Substituting in (13) we 
have 


dQ 

dt 


C-R 


Q 


(14) 


which has the same mathematical form as equations (2) 
and (12). As you can see in (14), the decay constant A, is 


given by 1/(CR) and must have units of seconds- 1 , jg 
Therefore, the half-life ln(2 )fk of the capacitor discharge is | 
the product of (C-R) and In(2). If you use a 200 ohm | 
resistor and a 5 volt, 1 farad capacitor, C-R = 200 seconds, j 
ln(2) = 0 .693, and the half-life is 139 seconds. Longer or (; 
shorter half-lives can be selected by changing the resistor I 
or capacitor. The capacitor can be charged with two 1.5 I 
volt batteries in series. Charge Q of the capacitor is I 
monitored during decay by measuring the voltage across I 
it with a multimeter. Even with the purchase of a I 
multimeter, the whole setup costs only about $60 - or you I 
can probably find all the parts in your local physics t 
department. The results are so good that in a plot of I 
voltage vs. time, the data are indistinguishable from an I 
exponential curve. 1: 

ALGEBRAIC EQUATIONS FOR 
RADIOACTIVE DECAY 

Students who understand algebra, but have not had m 
calculus, need not be disadvantaged by the mathematics ll 
of exponential decay. In addition to the analog fi 
experiments described in this paper, for many years I have 11 
used an algebraic derivation of the equations used for || 
radioactive dating. Even students who have had calculus |g 
tell me that they find this approach helpful. The -1. 
derivation is based on a simple recursive formula that is | 
developed as a discrete dynamical system, and is well- t 
suited to a spreadsheet. The recursion strategy used here fi 
is similar to many published explanations of compound 
interest, and the comparison of compounding annually, 
monthly, daily, etc. (e.g. Pierce, 2007; Moneychimp, 2008). 

A similar derivation for radioactive decay in an 11 
abbreviated form can be found in Guenther (1958). The f- 
following text is a formal presentation of the algebraic , 
equations, which 1 use as the basis for a class handout. 

We begin with the experimental observation that the 
number of radioactive atoms of one isotope that decay in f 
one unit of time is directly proportional to the number of 1 
radioactive atoms present. This means that the fraction I 
(k) of radioactive atoms of one isotope remaining after the I 
passage of one unit of time of length tk (e.g. one hundred § 
years) will be a unique value (k) that depends only on the I 
isotope (such as 40 K) and the length of the unit of time (tk). I 
In our analysis, this observation is expressed as a f) 
recursive formula, generalized, and then solved for the | 
total time (t) that has passed since the start of the decay I) 
process. 

Let Pu be the initial number of radioactive atom ("parents") i: 

in a sample. 1 

Let k be the fraction of radioactive atoms remaining after ft 
one unit of time. I 

Let tk be the length of the one unit of time. 

Let P n be the number of radioactive atoms ("parents") f 

remaining in the sample after n units of time, where n 
is an integer. jp 

After one unit of time, the number (Pi) of parents 1) 

remaining is: 1 : 

Pi = Pok (15) I 
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After two units of time, the number (P 2 ) of parents 

remaining is: 
o 

P 2 =Prk (16) 

Substituting for Pi from the first equation (15) we have: 

P 2 = (Po'k)-k = Po'k 2 (17) 


where X is a positive number that has units of (vr) -1 . The 
minus sign is used because k is a fraction and, therefore, 
ln(k) is negative. Using (26) to substitute for the term in 
brackets in (25) and changing notation (Pt = P„, where Pt is 
the number parent atoms remaining at time t) we have: 

{ \ 


In 



\Piii 


t-W 


(27) 


After three units of time, the number (P 3 ) of parents 


P 3 = P 2 -k (18) 

which becomes upon substitution for P 2 from equation 
(17): 

Pi = (Pirk 2 )-k = Pivk 3 (19) 

By continuing this procedure, we can show that after n 
units of time (tk) the number of parents remaining is: 

Pn = Pt)-k n (20) 

This equation expresses the fact that at anv time the 
number of parents remaining is related in a simple way to 
the initial number of parent atoms. For example, if we 
select our unit of time tk so that k = (1/2), then one-half of 
the parents will remain after one unit of time. This unit of 
time, ti/ 2 , which is different for each radioactive isotope, is 
called the half-life of an isotope. Using (20), we can see 
that after 4 half-lives (i.e. 4 units of time), the number of 
parents remaining is: 

1 \4 / 

1_ 

J6, 

As geologists, we are interested in using these 
relationships to find the total time (t) since the start of a 
decay process. To do so we start with a simple 
rearrangement of equation (20): 


I', • I’ 


V 4 


= Pn 


( 21 ) 


l Po > 


= k 


( 22 ) 


Taking the logarithm of both sides of the equation and 
using the identity In (a b ) = b In a we have: 


In 


VPo ) 


= if- In (k) 


(23 


The total time t (yr) is related to the number of units of 
time n and the length of one unit of time tk (yr) as follows: 

t = n ■ t k (24) 

Substituting for n in (23) using the relationship (24) and 
rearranging: 


In 


Pj_ 

fPo 




ln(k) = t 


Mk) 

t. 


(25) 


It has been shown experimentally that, for any one 
radioactive isotope, the term in brackets on the right of 
(25) has the same value whatever unit of time (t k ) is 
selected. This value is used to define a "decay 
constant" (1): 


X = - 


ln(k) 


U 


(26) 


Rearranging (27) to solve for t, we obtain: 


t = --In 

( \ 

Pl 

= — In 

1 \ 

Po 

X 

1P(U 

X 

ipj 


(28) 


which is the result we have been seeking. In these 
expressions, Po is the number parent atoms at time t=0, P t 
is the number parent atoms remaining at time t (vr), and X 
(\ r : ) is the decay constant. 

Our derivation of (28) assumes that t/t k is an integer n. 
This will not be correct for most times (t) if the value of k 
is a large fraction (e.g. V 2 ) or randomly chosen. However, 
because the unit of time tk is of arbitrary size, we can 
always choose the value tk so that i/tk is an integer. 
Therefore, equation (28) is true for all cases and, indeed, is 
the same as the exact solution obtained using calculus to 
solve the differential equation 


dP 

dt 


-A.P 


(29) 


that expresses, for a continuous dynamic system, the 
experimental observation we used to obtain (15). 


DISCUSSION 

It has been my observation that students remember 
more about their hands-on experiences in classes, labs, 
and field trips than about my lectures to them, and this 
trait has apparently- been true for Smith College students 
for a long time (Blakeslee, 1945). In this paper, I have 
compiled information about a number of hands-on 
laboratory experiments that I have found to be useful in 
teaching about exponential decay and radiometric dating. 
Any one of these experiments might be used if the 
teaching goal is simply to provide a physical analog to 
help students understand the mathematics of radiometric 
dating. Using several of the experiments together has the 
added benefit of making connections for your students 
among a number of physical processes they should know 
about as geologists. 

Throughout this paper I have emphasized the 
practical details of the decay experiments, providing 
advice on design and expected results (see Table 1 for 
observed half-lives). The learning that will result from 
adding these hands-on activities to your course depends 
on the classroom culture you create and the questions you 
ask. Don't give the students too much information; they 
will learn less if they are simply following your 
instructions. If you can divide the class into small groups, 
use some experiments that follow the exponential decay 
equation and others that do not. Challenge your students 
to prove the mathematical nature of the decay process 
they observe with graphical analysis. Have student 
groups trade experiments and apply or test each other's 
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TABLE 1. OBSERVED HALF-LIVES OF DECAY 
EXPERIMENTS 


Experiment 

Equipment 

Observed 

Half-Life 

(min) 

U-Tube Water 
Flow 

Two Plexiglass tubes, 1.75-inch i.d. 
Copper connecting tube, 0.5-inch 
i.d. 

0.3 - 67.4 

Siphon Water 
Flow 

Plastic tubing, 2 mm i.d., 4 feet long 
Connecting two vertical tubes, 1.75- 
inch i.d. 

13 

Cylindrical 
Water Flow 

Tennis ball tube, 6.5-inch high, 

3/32-inch hole 

1.3 


Tennis ball tube, 6.5-inch high, pin 
hole 

17 

Coffee Cooling 

Ceramic cup, 12-ounce, uncovered 

25 


Paper cup, 12-ounce, uncovered 

28 


Ceramic cup, 12-ounce, covered 

54 


Paper cup, 12-ounce, covered 

54 

Metal Cylinder 
Cooling 

A1 cylinder, 1.5-inch o.d., 3-inch 
long 

20 

Capacitor 

Discharge 

0.47 farad, across 335 ohm resistor 

1 farad, across 500 ohm resistor 

2 

6 


interpretation. Consider adding a mystery component 
such as finding the starting time of an already running 
experiment (e.g. Wise, 1990; AAAS Science Netlinks, 
2001b; Sunderman, 2007). When did the vessel start 
draining? When will the vessel be empty? When will the 
coffee be cool enough to drink? When was the cup of 
coffee poured? What information or assumptions are 
needed to answer these questions? Ask students to 
explain deviations from their predictions and then design 
tests of their explanations. Insist that uncertainties in data 
collection and manipulation are evaluated and 
considered. Give the students enough time and 
motivation to do a good job. Have fun! 

The skills students need for these experiments are 
ones that every scientist needs, including problem 
definition, thoughtful planning, careful lab work, good 
record keeping, appropriate mathematical analysis, critical 
thinking and interpretation, clear and complete written 
presentation of results. I believe this group of experiments 
is especially good for helping students reinforce or build 
their quantitative skills. The fact that a similar 
mathematical analysis can be applied to the many 
different processes emphasizes for students the utility of 
understanding exponential decay - and growth. Perhaps 
that knowledge will help them meet the challenges they 
face in a world where exponential growth of population 
and of resource consumption continue, in spite of good 
advice to the contrary (Bartlett, 1978). It may also help 


them appreciate logistic (or sigmoid) functions, which 
have been used by Hubbard and others to understand the 
consumption of finite resources and concepts such as 
"peak oil" (Hubbard, 1971; Deffeyes, K.S., 2003). 

The flows of water, heat, electricity, and atoms (by 
diffusion) are processes that can be described by similar 
differential equations. If the physical conditions are 
appropriate, the solutions to those equations are similar 
exponential decay equations like (28). This feature of 
these flows makes the experiments discussed here 
possible and, therefore, provides a possible starting point 
for more detailed consideration of those processes. In 
particular, if the boundary conditions do not represent 
flow between homogenous regions, other approaches are 
necessary to get solutions. The concepts and results of the 
analog experiments can be used to discuss and develop 
numerical, finite difference solutions for more 
complicated physical situations. Next steps could be 
simply mathematical, or they could be more experiments 
that use many vessels (e.g. Gilbert, 1979; Blanck and 
Gonnella, 2005) or circuits (Lawrence, 1970; Wunderlich 
and Peastrel, 1978). Alternatively, they could be to 
explore the underlying atomistic basis for the observed 
physics and the associated mathematics of probability and 
stochastic systems. 
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